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Hierarchical organization is a cornerstone of complexity and multifractality constitutes its central 
quantifying concept. For model uniform cascades the corresponding singularity spectra are sym¬ 
metric while those extracted from empirical data are often asymmetric. Using the selected time 
series representing such diverse phenomena like price changes and inter-transaction times in the 
financial markets, sentence length variability in the narrative texts, Missouri River discharge and 
Sunspot Number variability as examples, we show that the resulting singularity spectra appear 
strongly asymmetric, more often left-sided but in some cases also right-sided. We present a unified 
view on the origin of such effects and indicate that they may be crucially informative for identifying 
composition of the time series. One particularly intriguing case of this later kind of asymmetry is 
detected in the daily reported Sunspot Number variability. This signals that either the commonly 
used famous Wolf formula distorts the real dynamics in expressing the largest Sunspot Numbers or, 
if not, that their dynamics is governed by a somewhat different mechanism. 
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Multi-scale approach aims at bridging the wide 
range of time and length scales that are inherent in a 
number of essential processes in complex natural phe¬ 
nomena. Multifractality constitutes a leading concept 
towards quantifying the related characteristics i- By 
now it finds applications in essentially all areas of the 
scientific activity including physics biology fij-lil , 

chemistry [10,[lT|, geophysics [1 21 113). e conomics |14| - l2l| . 
hydrology j , atmospheric physics [23]], quantitative lin¬ 
guistics miHj] , behavioral sciences [26j. music [27l l28j. 
and even ecological sciences [29|. 

At present the most efficient, numerically stable 
and precise [30] method to quantify multifractality is 
based on the Multifractal Detrended Fluctuation Anal¬ 
ysis (MFDFA) HU. Accordingly, for a discrete signal 
N one starts with the signal profile X(j) = 

Yjl = i < x >)> j = 1) N, where < ... > denotes 

averaging over all i’s. Then one divides the A'(j') into 
M s non-overlapping segments of length s (s < N) start¬ 
ing from both the beginning and the end of the signal 
(there are thus in total 2 M s segments). For each seg¬ 
ment a local trend is estimated by fitting an Zth order 
polynomial P„\ which is then subtracted from the sig¬ 
nal profile. For the so-detrended signal a local variance 
F 2 ( v , s) in each segment v is calculated for the scale vari¬ 
able s. Finally, by averaging F 2 (y , s) over all segments v 
one calculates the qth order fluctuation function: 

( 2M S . l/q 
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and q G R. The optimal range is q G [—4,4] [32]. 

The scaling behavior of F q (s) ~ s indicates fractal 
structure with the singularity spectrum f(a) jl] 

f(a) = q[a - h(q)] + 1 (2) 

for a = h{q ) + qh'(q). If h{q) = const, the signal is 
monofractal. A nontrivial g-dependence of h(q) indicates 
its more convoluted fractal organization termed multi¬ 
fractal. For model multifractal series /(a) typically as¬ 
sumes shape of an inverted symmetric parabola while the 
empirical ones often develop asymmetries. Here we ex¬ 
plore this issue and document that asymmetries in /(a) 
may provide important information about the series or¬ 
ganization. 

Spectacular examples, largely novel, of unques¬ 
tioned asymmetric multifractal spectra, generated using 
MFDFA algorithm, for signals representing diverse, mu¬ 
tually remote areas are shown in Fig. [1] They are seen 
to even display two distinct kinds of asymmetry. The 
DAX and DJIA stock market indices, the Forex market 
represented by the three exchange rate pairs (CHF/JPY, 
GBP/CHF and GBP/JPY), the series representing the 
sentence length variability in some selected literary texts 
(here Rayuela by Julio Cortazar and 2666 by Roberto 
Bolano; most of the world famous literary works are es¬ 
sentially monofractal in the sentence length variability 
and just a dozen or so are convincingly multifractal and 
the two considered here at the same time belong to those 
few that develop the most asymmetric /(a)) and the dis¬ 
charge of Missouri River [§3} show pronounced left-sided 
asymmetry. A rarer, right-sided asymmetry is displayed 
by the other two examples, the stock market company 
(here ADS and BMW from DAX) inter-transaction times 
and the Sunspot Number variability [H|. In order to es¬ 
timate a possible additional contribution to asymmetry 
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FIG. 1. (Color online) Singularity spectra f(a) for (a) daily returns of the two (DJIA and DAX) stock market indices over 
the period January 12, 1990 - October 12, 2013 (5881 data points), (b) minutely Forex returns on selected pairs (CHF/JPY, 
GBP/CHF, GBP/JPY) of currencies over the period 9:00 pm January 2, 2004 - 9:00 pm 20 March, 2008 (1703520 data points 
each), (c) the series representing length, expressed in terms of the number of words, of the consecutive sentences in the two 
narrative texts: Rayuela by Julio Cortazar (9848 sentences) and 2666 by Roberto Bolano (21319 sentences), (d) hourly changes 
of the Missouri River discharge (in cubic feet per second) at Waverly, taken from the U.S. Geological Survey database [33] over 
the period October 1, 1987 - October 1, 2007 (213344 data points), (e) inter-transaction intervals over the period November 28, 
1997 - December 31, 1999, for the two German stock market companies: ADS (167208 data points) and BMW (173664 data 
points) (f) daily Sunspot Number over thejperiod January 1, 1900 - June 30, 2014 (41817 data points); Source: WDC-SILSO, 
Royal Observatory of Belgium, Brussels [33]. The insets show the corresponding F q (s ) for q G [—4,4] with the white strip 
indicating q = 0. The detrending polynomial used, as optimal, is of the second order. In case (f) the upper limit of the scaling 
parameter s is less than half of the 11 years solar activity cycle which makes this procedure unbiased by the periodic trend. 


of the potential ’oscillating singularities’ (35| we also ap¬ 
plied the so-called wavelet leaders algorithm [36|, |3^] to 
all these data but detected no significant signals of such 
oscillations. 

The left side of f{a) is determined by the positive q- 
values, which filter out larger events and the opposite 
applies to its right side. Hence, asymmetry in f(a) sig¬ 
nals non-uniformity of the underlying cascade. The (a)- 
(d) cases in Fig. |T] are thus seen to be more multifrac¬ 
tal in arrangement of the large events and far less such 
in the small ones. In some of those cases (CHF/JPY, 
GBP/CHF or 2666) the right-side of /(a) contracts so 
strongly that it indicates essentially a monofractal char¬ 
acter of the corresponding small fluctuations. This sug¬ 
gests that the entire signal can be considered a mix¬ 
ture of the large scale homogeneous multiplicative cas¬ 
cade, giving thus rise to a pronounced left wing in /(a), 
and of the small scale noise-like background, shrinking 
the corresponding right wing in /(a). Such a noise-like 


component in empirical data is not very unusual as it 
may for instance originate from the measurement uncer¬ 
tainty or from some coarse graining that effectively af¬ 
fects more just the small scale fluctuations. The reverse 
kind of asymmetry - right-sided - applies, however, to 
the cases (e) and (f), and here the situation is more in¬ 
triguing and unusual since it indicates that pronounced 
multifractality operates on small scale fluctuations while 
the dynamics of large scale fluctuations is much poorer 
in this respect. That the stock market company inter¬ 
transaction intervals may be governed by such a dynam¬ 
ics is conceivable. Small intervals occur between trans¬ 
actions grouped within the same clusters of an enhanced 
volatility on this particular company and may thus be 
strongly, also nonlinearly correlated. The distances be¬ 
tween such in time more separated clusters of activity, 
thus larger inter-transaction intervals connecting them, 
may be less correlated [32j |. 

More difficult to interpret - thus even more interesting 
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FIG. 2. (Colore online) Examples of /(a) for the two binomial 
cascades with A = 1.05 and A = 1.5, the average / a „(a) of 
such two /(a)’s and /„(a) calculated for the sum of these 
two cascades. The upper two insets show explicit spreads of 
these two cascades. To better reflect proportions in the same 
scale the original values of A = 1.05 cascade are multiplied by 
a factor of 50. The lowest inset displays F q (s) functions for 
the sum of both cascades. 


- is the analogous result for the Sunspot Number vari¬ 
ability. The smaller fluctuations (seen through q < 0) 
develop a broad right-sided f(a) which suggests their 
cascade-like hierarchical organization. At the same time 
the large ones, filtered out by q > 0, show at most a rem¬ 
nant of multifractality. They thus do not belong to the 
same hierarchy or, at best, strongly distort it. This ob¬ 
servation may either demand revision of the famous Wolf 
formula (R = k* (N s + 10 N g ), N s - number of spots, N g - 
number of groups) 1.38I | commonly used for expressing the 
Sun activity in terms of the number of spots or may indi¬ 
cate that the corresponding large and small activities are 
governed by a somewhat different cascading mechanisms. 

The quality of scaling can be assessed from the insets to 
the corresponding panels of Fig. Q] which show the F q (s) 
dependencies on s. In (a) - (d) one in addition sees that 
the scaling exponents (slope of F q (s )) visibly depend on 
q for q > 0 while very weakly for q < 0. The opposite is 
true for (e) and (f) cases. Naturally, this correlates with 
orientations of asymmetries seen in f(a). 

The most likely origin and thus interpretation of the 
above effects can transparently be illustrated by using 
the model cascades. Here we use the binomial cascade 
with multipliers whose logarithms are drawn from the 
Gaussian N( A, 1) distribution. Examples of f(a) for two 
such cascades with A = 1.05 and A = 1.5, the average 
/ a „(a) of such two /(q')’s and f SU m{ot ) calculated for the 
sum of these two cascades are shown in the main panel of 
Fig. [2] While / at ,(a) is symmetric, the f sum (ct) is much 
narrower and already strongly asymmetric. To quantify 
these characteristics we introduce the asymmetry param- 
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FIG. 3. (Color online) Singularity spectra of /(a) for an in¬ 
creasing number of the superimposed binomial cascades gen¬ 
erated independently but all with the same A = 1.1. Inset 
illustrates how the corresponding asymmetry parameter A a 
evolves. 


eter 


A a = (A a L - A a R )/{Aa L + Aa R ) (3) 

where Aa L = a 0 - a min and A a R = a max - a 0 and 
oimin, ctmax , OiQ denote the beginning and the end of f(a) 
support, and the a value at maximum of f(a) (which cor¬ 
responds to q = 0), respectively. A simple superposition 
of only two such stochastic binomial cascades with differ¬ 
ent parameters produces strongly left-sided asymmetric 
(A a = 0.55) multifractal spectrum /(a), resembling the 
corresponding empirical ones from Fig. [T] It is instruc¬ 
tive to see that the left side of f SU m(ct ) coincides with 
f(a) for the A = 1.5 cascade while its right side with 
the one for A = 1.05 cascade. The origin of this result 
becomes clear from the relative proportions of the size 
of fluctuations in these two individual cascades as seen 
in the upper two panels in Fig. [2] The A = 1.05 fluc¬ 
tuations are on average about two orders of magnitude 
smaller than the A = 1.5 ones. Thus in their sum the 
small fluctuations are dominated by A = 1.05 and fil¬ 
tered out by q < 0 which corresponds to the right side in 
fsum and the opposite applies to large fluctuations. An 
equivalent way to see this effect is through F q (s) defined 
by Eq.|T]and shown in the lowest panel of Fig. [2] Indeed, 
all F q (s) scale with s but spread much more for q > 0, 
thus width on the left side of f sum {a ) is larger. 

Similar effect of left-sided asymmetry one even obtains 
by superimposing cascades generated with the same A. 
The result for an increasing number N of independently 
drawn and then superimposed cascades with A = 1.1 is 
displayed in Fig. [3] Here however more components are 
needed to reach the same degree of asymmetry as in the 
Fig. [2] case. The left-sided asymmetry originates here 
from the fact that small fluctuations are more abundant 
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FIG. 4. (Color online) Singularity spectrum f(a) (i) for 
A = 1.05 binomial cascade and (ii) for the one obtained from 
it by setting a threshold of 4.7, then from all events that ex¬ 
ceed this threshold value subtracting the over-threshold value 
(which affects 1% of all events) and finally by adding consecu¬ 
tively to those affected one number drawn from a normal dis¬ 
tribution, one independently to each. Upper inset, from top 
to bottom, shows the original binomial cascade with the trun¬ 
cation threshold (dashed line), the white noise added and the 
resulting final series. Lower inset displays the g-dependence of 
the fluctuation functions F q ( s ) for the so-obtained final series. 


in such individual cascades therefore summing up an in¬ 
creasing number of them one sooner approaches the white 
noise limit, i.e., destroys the original hierarchical organi¬ 
zation on the level of small fluctuations than on the level 
of the large ones. The right side of the resulting f(a ) 
contracts sooner to a monofractal with an increasing N 
than does the corresponding left side. This can be seen 
from Fig. [3l especially from the inset which shows the 
A a dependence on N. After the initial sharp increase it 
reaches maximum at around N = 30 with A a « 0.7 and 
then slowly starts decreasing. For N = 1000 the f(a) 
already is much closer to a monofractal with A a ss 0.2. 

While the left-sided multifractal asymmetry easily 
emerges when superimposing cascades, the right-sided 
one is more peculiar to model. From an algorithmic per¬ 
spective it implies just a reverse,i.e., a more uniform hier¬ 
archical organization on the level of smaller fluctuations 
and more noise-like behavior of the large fluctuations. 
Such an interpretation indicates construction of a series 
obeying condition of this kind. An example is displayed 
in Fig. [31 This is the binomial cascade generated from 
the same A = 1.05 model as before, but now a top sec¬ 
tion of the largest events is randomized as follows. 1% 
of the largest events is selected which in this case cor¬ 
responds to those that exceed a threshold of T = 4.7cr. 
Then, their values are replaced by the sum of T and of a 
number drawn independently for each event from 7V(0,1). 


The net result of such a construction preserves organiza¬ 
tion of the original small fluctuations and randomizes the 
largest ones. The so obtained series is displayed in the 
upper inset of Fig. [3] (final series), together with the two 
components that are used to form it. The correspond¬ 
ing F q (s) functions (lower panel in Fig. |4j) become now 
more spread for q < 0, thus the f(a) spectrum becomes 
right-sided (A a = —0.4), which confirms the demanded 
reversed organization. Of course, by playing with T or 
with noise added one may modulate the degree of this 
right-sided asymmetry. 

In conclusion the effects of asymmetry in singularity 
spectra f(a) of the time series representing complex sys¬ 
tems may not just be numerical artifacts as they seem to 
be treated in the literature on the subject, but instead 
may contain genuine information about the composition 
of such series.They then indicate at least two different 
regimes in the multifractal scaling. The left-sided asym¬ 
metry corresponds to a pronounced multifractality on the 
level of larger fluctuations and its suppression towards 
monofractality when going to small fluctuations. These 
small ones thus resemble a noisy background accompany¬ 
ing the entire signal. Globally, such structures can easily 
be modeled as a superposition of independent, multifrac- 
tally symmetric cascades. Even in relation to reality such 
a modeling seems appropriate since participation of sev¬ 
eral somewhat independent multiplicative factors on a 
given phenomenon seems indeed possible. Out of those 
studied here the stock market index is an explicit sum 
of prices of the individual companies entering the basket. 
The span of f(a) of such sums appears also visibly nar¬ 
rower. This thus calls for caution when interpreting the 
width of f(a) as a measure of the ’degree of complexity’. 
In case of asymmetric f(a) application of an appropriate 
multifractal cross-correlation analysis [40l - l42| in order to 
disentangle the original signal into uniform cascade com¬ 
ponents may in this connection be recommended. In the 
right-sided asymmetry the situation is reversed. Here, 
these are the smaller fluctuations that develop a pro¬ 
nounced multifractal hierarchy while the largest ones are 
suppressed. In relation to natural phenomena this seems 
to be a much rarer category of distortion. Especially in¬ 
triguing is the Sunspot Number variability case belonging 
to this category. The arguments presented here indicate 
that possibly the famous Wolf formula 38] commonly 
used to digitally express this aspect of the Sun activity 
does not represent a faithful mapping for the large events 
and needs to be refined such that the character of cas¬ 
cading from the smaller events be extended to the large 
ones. Of course, at the present stage another possibility 
that the large events are governed by a somewhat dif¬ 
ferent dynamics cannot be excluded either. In any case, 
resolving this issue emerges as a great scientific challenge. 

We thank Drs Frederic Clette and Laure Lefevre of the 
Royal Observatory of Belgium, Brussels, for very helpful 
exchanges on the issue of the Sunspot Number variability. 
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